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ABSTRACT 


The investigation dealt with the numerical and experimental 
investigation of the effect of porosity on the vortex shedding from a 
cambered plate. The effect of porosity was incorporated into the anal- 
ysis through use of the experimentally obtained pressure-drop data. 
The results have shown that the porosity of the existing materials used 
in parachute canopies is not sufficient to eliminate the negative pres- 
sure gradient. Canopies with nonuniform porosities will have to be 
designed to overcome the collapse phenomenon resulting from the 


unfavorable pressure gradients. 
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I. INTRODUCTION 


The current parachute-canopy design methodology is mainly 
empirical. It is based on either model tests in wind tunnels (see, e.g., 
Klimas and Rogers, 1977; Pepper, 1986; Rychnovsky, 1977) or on the 
deployment of full-scale parachutes with test loads. However, under 
certain conditions, parachutes have been observed to collapse during 
periods of rapid deceleration, subsequent to the initial inflation of the 
canopy with a resulting loss of drag and ultimately the payload. 

Studies began with a two-dimensional model of flow about a rigid, 
nonporous camber. Even though highly idealized, the purpose of the 
investigation was not a solution of an immediate practical problem but 
to understand the physics of the problem. With that in mind, the dis- 
crete vortex model (DVM) (see, e.g., Chorin, 1973; Sarpkaya, 1975) 
was used to model the flow past bodies with included angles of 120, 
180, and 240 degrees, and was analyzed by Mostafa (1987) and Munz 
(1987). These investigations have shown that the drag force is pre- 
dicted fairly well during the steady-state period and up to the middle 
of the deceleration period. Beyond that time, the predicted forces 
tended to be relatively larger than those measured experimentally. 

Nevertheless, the forces measured during the final stages of 
deceleration were very small, so the degree of exactness was not a 
concern. Moreover, the reasons leading to the collapse (backwash of 


the vortices) occurred shortly after the start of deceleration, where 


the predicted and measured forces compared well. Previous 
investigations at the Naval Postgraduate School have concluded that 
the development of negative differential pressures in the central 
region of the camber is the likely cause of parachute collapse at high 
deceleration rates (see, e.g., Munz, 1987). 

In view of the foregoing, it became clear that the flow through the 
porous camber, if managed properly, can help to minimize or elimi- 
nate the collapse phenomenon by eliminating the negative differential 
pressure across certain sections of the camber. The present 
investigation was undertaken to study in detail the effect of porosity on 
a two-dimensional rigid camber in order to bring the theoretical 
problem a step closer to reality. In the following, the fundamental 
concepts used in the introduction of flow through a porous wall and 
the results obtained with various methods and porosities are described 


in detail. 


Il. EFFECTS OF POROSITY 


A. REVIEW OF THE FLOW THROUGH POROUS SCREENS 

A parachute canopy is a porous material. The amount of flow that 
goes through such a material depends on the percentage area opening, 
the shape and size of the threads, the moisture content of the mate- 
rial, the differential pressure across the canopy, the angle of approach 
of the ambient flow, the magnitude and direction of the flow on the 
downstream side of the canopy, and the degree of time-dependence of 
the upstream and downstream flows (i.e., the vortex motion down- 
stream of a parachute). 

In order to incorporate the effect of porosity into the analysis, one 
had to know at least the flow rate through the canopy as a function of 
differential pressure, assuming that all of the other parameters cited 
above were of secondary importance, or that their effects may be 
incorporated into a flow coefficient relating differential pressure to 
the normal velocity of the flow rate. This seemingly simple search for 
a reliable (theoretical and/or experimental) equation turned out to be 
anything but simple. The following is a brief description of the previ- 
ous investigations of the flow past screens. 

A fairly detailed review of the flow through screens is given by 
Laws and Livesy (1978). There are a number of investigations (see, 
e.g., Baines and Peterson, 1951; Turner, 1969; Reynolds, 1969; 
Graham, 1976) which dealt with the flow past fairly porous screens for 


the purpose of understanding, managing, or generating isotropic tur- 
bulence. Often the resistance through the screen was of minor 
importance. Others (see, e.g., Taylor, 1944; Taylor and Davies, 1944; 
Wieghardt, 1953) who investigated the flow through infinite screens 


expressed the resistance as 
Ap = 7pKU" 


Koo and James (1973) and Cumberbatch (1982) investigated the flow 
past finite screens set in infinite flows, but the effects of body shape, 
Separation, vortex formation, flow magnitude and direction down- 
stream of the body, and deceleration of the flow were not studied. It 
became clear that there was an enormous amount of research done on 
the flow through screens and yet there was very little information 
which could be extracted from them for use in this study of a deceler- 
ating parachute. Payne (1978) reviewed the data from many 
researchers in an effort to relate the geometric porosity (open area 
ratio) to the differential pressure, and developed an equation in terms 


of dynamic pressure loss term and a viscous loss term, 
Ap = K,U* +K,U 


However, this study involved many different cloth types, not just 
parachute cloth, and included no evaluations of the effects of flow 
direction. This prompted the undertaking of an experimental investi- 


gation to establish the proper relationship between the pressure drop, 


—_——~ 





the flow rate, and the flow direction on the actual samples of canopy 


materials used in parachutes. 


B. PARACHUTE CLOTH EXPERIMENTS 
1. Experimental Appar nd Procedur 

The purpose of the experiments was twofold: to find the dif- 
ferential pressure versus flow rate relationship, and to find the 
change, if any, of this relationship with the downstream flow direction. 
For this purpose, an experimental apparatus was conceived, designed, 
and manufactured as shown in Figures 1 and 2. The flow rates at 
points A and B were measured with the standard flow meters, and the 
static pressures upstream of the screen were recorded for each flow 
rate through the screen. Static pressure downstream of the screen 
was assumed to be atmospheric. 

Observations have shown that the canopy sample did not 
bulge measurably during the experiments, even at the differential 
pressures encountered by a parachute beginning to descend at a speed 
of about 600 miles per hour. This was primarily due to the proper 
construction and sandwiching of the circumference of the canopy 
sample. However, this is not meant to exclude the possibility that at 
high pressures, the slight cambering of the sample may have stretched 
it and thereby increased its porosity slightly, leading to a smaller dif- 
ferential pressure for a given flow rate than that had it not been 
Stretched. In fact, this was checked by conducting the experiments 


first with increasing flow rate and then decreasing the flow rate. No 





Figure 1. Parachute Cloth Differential Pressure 
Measurement Apparatus, Front View 
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Figure 2. Parachute Cloth Differential Pressure 
Measurement Apparatus, Three-Quarters View 


attempt was made to wait a long time after the completion of the 
experiment with ascending pressures to see if the material would 
return to its original unstretched shape, or whether fibers used in 
construction of the canopies develop a permanent set beyond certain 
differential pressures and retain their deformed shape. This leads to 
the possibility that a once-used parachute may fall faster than a never- 
used one. 

Of the seven available canopy materials, four samples were 
chosen and tested to provide a range of data. The characteristics of 
the four samples are shown in Table 1. The identifying letters are 


used for later reference. 


TABLE 1 
CHARACTERISTICS OF PARACHUTE CLOTH SAMPLES 


Military 
Sample Specifications 
Letter Material Weight Number 
A Nylon 1-1 oa = 
B Nylon 4.75 Oz. MIL-C-8021D 
Type | 
C Nylon 7.0 Oz. MIL-C-8021D 
Type Il 
D Kevlar SEO) (er4, _ 


During the course of the experiments, the flow on the down- 


stream face of the sample was systematically changed to see if the 


differential pressure across the sample would increase for a given flow 
rate through the sample. 
2. Results 

The results of tests on the effect of flow on the downstream 
face of the samples have emphatically and somewhat surprisingly 
shown that the flow through the screens investigated is not affected by 
the presence or absence of the downstream flow. This apparently 
simple finding is of extreme mathematical and practical importance. 
It simplifies the analysis and the possible management of porosity to 
eliminate parachute collapse. 

The results obtained with the four representative canopy 
samples are shown in Figures 3 and 4. 

Figure 3 shows the variation of differential pressure as a func- 
tion of the dynamic head, tpV. , where Vn is the velocity through the 
canopy (2). It is clear that the value of differential pressure varies 
with the weight of the sample and manufacturing. For samples B and 
C, the military specifications expressed the mean velocity as flow rate 
per unit area. For both samples, the military specifications call for 50 
to 90 cfm/ft2 at a differential pressure of .5 inch of water, and 450 to 
650 cfm/ft2 at a differential pressure of 20 inches of water. Sample B 
extrapolated to 88 cfm/ft2 and 654 cfm/ft2, respectively, and sample C 
extrapolated to 89 cfm/ft2 and 487 cfm/ft2, respectively. 


Figure 4 shows the pressure coefficient Cp, defined by 


2 Ap 
C,= —, 
P pv 
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Figure 3. Differential Pressure vs. Dynamic Head 
for Parachute Cloth Samples 
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Figure 4. Pressure Coefficient vs. Dynamic Head 
for Parachute Cloth Samples 
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for all the samples as a function of dynamic head. Clearly, the pressure 
coefficient remains nearly independent of velocity for a given canopy. 
As noted earlier, the flow on the downstream face of the sample 
(which ranged from O to 14.4 ft/sec) had no effect on the results 
Shown in Figures 3 and 4. Table 2 shows the minimum, maximum, 
and average values of the pressure coefficient for the different 


samples. 


TABLE 2 


COMPARISON OF:-THE PRESSURE COEFFICIENTS FOR 
THE DIFFERENT CANOPY MATERIAL SAMPLES 


Sample Minimum Cp Maximum Cp Average Cp 
A 2003 27.3 (32.4) 22.9) (24s) 
B 518) 10; 43.0 (52.0) 36.7 (39.2) 
C 39.0 48.9 (54.8) 42.8 (44.8) 
D 61,5 73.3 (80.4) 66.3 (69.1) 

Notes: 


(1) The values in parentheses in the maximum column 
are for the lowest velocity data points. 


(2) The average values do not include the pressure coeffi- 
cients for the lowest velocity data points because they 
were deemed not representative of the velocities that 
parachutes descend. The values in parentheses in the 
average column are averages if these data points had 
been included. 
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Hit. ER PE ENT 


A. THEORETICAL FOUNDATIONS OF THE MODEL 

The flow about a cambered plate is obtained by successive trans- 
formations from a circle plane. In fact, it can be shown that the flow 
in the C-plane (circle plane) may be transformed to the flow in the 
z-plane (camber plane) through the use of 


b= 2m?-1 


DN ae 71 


in which m is a parameter related to the total angle of the circular arc. 
For example, m equal to .5, .707, or .866 corresponds to camber 
angles of 120, 180, or 240 degrees, respectively. Figure 5 shows the 


transformation from a circle to a 120-degree camber. 


, & 
i 60 
| 


Figure 5. Circle and Physical Planes 
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The complex function describing the flow is given by 








2 it, ily, 2 
w=-u[t+$]+ L(G = coe Ln -t73 
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in which yp and Cxp represent respectively the strength and location 
of the k-th p-vortex (counter-clockwise rotating), and xq and (xq 
represent the strength and location of the k-th q-vortex (clockwise 
rotating). As described in Munz, discrete vortices are introduced at 
the tips of the camber at regular time intervals, and evolution of the 
vortex wake is calculated through the use of the velocities induced at 
each vortex location. Furthermore, the forces acting on the camber 
are calculated through the use of the generalized Blasius theorem, or 
through the integration of the differential pressure. The numerical 
model requires the input of a velocity history, and in the present 
investigation, the velocity profile that was employed was that used by 
Munz for the simple purpose of delineating the effect of porosity 
relative to the non-porous case calculated by Munz. 

The existing analysis had to be considerably modified in order to 
incorporate the effect of porosity. The inclusion of the effect of 
porosity was already shown to be difficult as far as the previous 
experimental investigations are concerned. It turns out that its math- 


ematical analysis is just as difficult as its experimental understanding. 
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B. INOQUE’S METHOD 

The only work which dealt with a steady uniform flow past a sim- 
ple porous flat plate was carried out by Inoue. To emphasize, the flow 
was not subjected to any deceleration and the only unsteadiness intro- 
duced into the flow came from the evolution of the vortex shedding 
following the impulsive start of the flow at a constant velocity. Inoue 
made no attempt to relate the differential pressure across the plate to 
the flow rate. According to his method, the complex function for a 


plate normal to the flow is given by 


o(6-3) a-0(6+4) 
2 


2 
N _ —Ir, §€-¢ 
+1Dphing > riXglne = 








in which a is an artificial porosity coefficient related to the geometric 
porosity and aUjnf is the part of the flow which passes through the 
plate uniformly, regardless of the differential pressure distribution 
prevailing along the plate. Inoue did not calculate forces or pressure 
distributions, and confined himself to the comparison of the calculated 
flow kinematics with sample flow pictures, presumably obtained under 
conditions similar to those encountered in the numerical calculations. 

Inoue’s analysis, though not applicable to the analysis of the actual 
flow under consideration, was tried anyway, to gain some experience 
regarding the overall performance of the code developed in the pres- 
ent investigation. Figures 6 through 10 show, at representative times, 


the evolution of flow for a = 0.25 and for the velocity profile given by 
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Figure 6. Wake Vortex Positions for a = 0.25 at T = 8.85 
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Figure 7. Wake Vortex Positions for a = 0.25 at T = 11.29 
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Figure 8. Wake Vortex Positions for a = 0.25 at T = 12.92 
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Figure 9. Wake Vortex Positions for a = 0.25 at T = 14.54 
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Figure 10. Wake Vortex Positions for a = 0.25 at T = 16.98 
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U7 1 for i Uot/c Ss 9.72 


and 


U/Up, = 0.97T - 0.05T2 - 3.70 


and 


A = 0.97 - 0.10T 


in the interval 9.72 < T < 11.48 and 


U/Uo = - 0.3423 + 0.0072T2 + 3.82 


and 


A = - 0.3423 + 0.01445T 


in the interval 11.48 < T < 17.95. 

For T greater than 18, the velocity and acceleration are zero. 

Figures 11 through 15 show the flow at exactly the same times 
about a nonporous camber (a = 0). Comparison of these figures shows 
that, as anticipated,the effect of the porosity is to push the wake fur- 
ther downstream and thereby minimize the consequences of the 
return of the wake to the body. In fact, Figures 16 through 20 show 
the comparison of the differential pressures at the corresponding 
times for a = 0.0 and 0.25. Apparently, the differential pressure for a 


= .25 becomes less negative than the one for a = 0.0, even though, as 
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Figure 11. Wake Vortex Positions for 
Non-Porous Camber at T = 8.85 
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Figure 12. Wake Vortex Positions for 
Non-Porous Camber at T = 11.29 
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Figure 13. Wake Vortex Positions for 
Non-Porous Camber at T = 12.92 
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Figure 14. Wake Vortex Positions for 
Non-Porous Camber at T = 14.54 
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Figure 15. Wake Vortex Positions for 
Non-Porous Camber at T = 16.98 
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Figure 16. Differential Pressure Comparison at T = 8.85 


ZH, 


DP 
1 
Se he 


60. #501749" CAMGe eee sete 0 10 20 30 40 50 60 


Nonporous Camber 





Porous Camber, a = 0.25 


Figure 17. Differential Pressure Comparison at T = 11.29 
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Figure 18. Differential Pressure Comparison at T = 12.92 
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Figure 19. Differential Pressure Comparison at T = 14.52 
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Figure 20. Differential Pressure Comparison at T = 16.98 
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noted earlier, Inoue’s method does not take into account the 
relationship between the flow and the differential pressure, or the 
variability of a along the camber and with time. In other words, a 
must vary with time and with the angle along the camber. Thus, the 
comparison of the a = 0.0 and a = 0.25 cases will not be made in 


greater detail than Inoue’s method deserves. 


C. PRESENT METHOD OF POROSITY MODELING 

With the realization of the fact that the flow through the camber 
must be related to the prevailing differential pressure and that the 
direction of the flow in the wake has no measurable effect on the dif- 
ferential pressure versus dynamic head, a new method had to be 
devised which will properly account for the effect of porosity. 

The present method takes advantage of the simplicity of Inoue’s 
method and accounts for the local dependence of the flow on 
differential pressure. For this purpose, the value of a is chosen so as 
to make the Vy at @ = 0.0 proportional to the / Ap as given by Figure 
4. The variations in the flow across the camber beyond and above that 
given by Inoue’s method are accounted for through the use of discrete 
vortices placed along the camber. A number of sample calculations 
have shown that the use of 118 body vortices (one every degree) on a 
120-degree camber will be more than sufficient to account for the 
additional cross-flow effects smoothly. The first vortices were placed 
one degree from the tips, and no vortex was placed at 8 = 0. The mid- 


points between the vortices were used as the control points to 
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calculate the velocity normal to the camber. At each control point, the 


differential pressure was calculated, and then Vy was found from 


_ 1 /24Ap 
Yo 


The velocity induced by the uniform flow based on a was subtracted, 
and then the strengths of the body vortices were determined by 
equating the remaining velocity to that induced by the body vortices at 
the control points. Once the strengths of the body vortices were 
determined, the velocity induced at each vortex in the wake was cal- 
culated, the velocity of the ambient flow was updated, and all of the 
wake vortices were convected for a time interval At using a 274 order 
modified Eulerian scheme. Then the procedure was repeated, up to 
any specific time (in this analysis T = 24). The wake vortices were 
never merged together, no matter how close they came, and were only 
removed from the simulation if they contacted the camber. To pre- 
vent abnormally high induced velocities from occurring due to the 
proximity of one vortex to another, whether body or wake vortices, 
Rosenhead’s (1930) smoothing scheme was used, where the complex 


velocity at a point was multiplied by 


x? 


x48 





where x is the distance from the point to the vortex inducing the 
velocity and 6 is the attenuation factor. In this analysis, 6 was chosen 


to be 0.05. 
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The method described above is quite general and could be used to 
investigate the effect of the porosity of any canopy. Because it is 
deemed sufficient to understand the physics of the phenomenon, only 
two canopy characteristics were used in the numerical experiments. 
The first case corresponded to the characteristics of canopy sample A 
(Cp = 22.9, see Table 2), and the second to a canopy which is much 
more porous than any of the samples tested, with an assumed Cp = 
4.45, which corresponded to an initial value of a = 0.25 at 6 = O and 
=>: 

Figures 21 through 25 show the evolution of the wake for Cp = 
22.9 (sample A). Comparison of these figures with Figures 11 through 
15 shows that the effect of porosity, however noticeable, is not large 
enough to show a significant downstream displacement of the wake 
vortices. The importance of this conclusion lies in the fact that the 
most porous parachute canopy used by the industry is not sufficiently 
porous to minimize the consequences of the wake return, i.e., the 
collapse of the parachute as evidenced by field experiments. Thus, it 
has been possible through a relatively idealized analysis to show the 
inadequacies of the materials used in the construction of large 
parachutes. 

The next step in the analysis was to show the beneficial effects of 
porosity when it can be made sufficiently large. For this purpose, a 
value of Cp = 4.45 was chosen, for the reasons noted earlier. Figures 


26 through 30 show the evolution of the wake and Figures 31 through 
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Figure 21. Wake Vortex Positions for Cp = 22.9 at T = 8.85 
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Figure 22. Wake Vortex Positions for Cp = 22.9 at T= 11.29 
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Figure 23. Wake Vortex Positions for Cp = 22.9 at T = 12.92 
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Figure 24. Wake Vortex Positions for Cp = 22.9 at T = 14.54 
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Figure 25. Wake Vortex Positions for Cp = 22.9 at T = 16.98 
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Figure 26. Wake Vortex Positions for Cp = 4.45 at T = 8.85 
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Figure 27. Wake Vortex Positions for Cp = 4.45 at T = 11.29 
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Figure 28. Wake Vortex Positions for Cp = 4.45 at T = 12.92 
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Figure 29. Wake Vortex Positions for Cp = 4.45 at T = 14.54 
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Figure 30. Wake Vortex Positions for Cp = 4.45 at T = 16.98 
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35 the differential pressures. It is evident from a comparison of Fig- 
ures 26 through 30 with either Figures 11 through 15 (nonporous 
case) or Figure 21 through 25 (Cp = 22.9 case) that the wake is 
significantly displaced downstream and that, as seen in Figures 31 
through 35, the differential pressure is far less negative than those 
shown in Figures 16 through 20 (a = .25 and a = 0.0). Figure 36 
shows histories of the differential pressure computations at 6 = 0 


22.9, 


degrees for the nonporous case, the present model with Cp 
Inoue’s method with a = 0.25, and the present model with Cp = 4.45. 
It is easily seen that the present method with Cp = 22.9 has differen- 
tial pressures only slightly less in magnitude than the nonporous case, 
while the two other, more porous cases have differential pressures 
significantly smaller in magnitude and delay the occurrence of nega- 
tive differential pressure at the center of the body. It is also seen that 
while the two large porosity cases during steady state and the early 
deceleration periods are very similar, during later stages of the decel- 
eration the present method results in lower differential pressures as it 
accounts for the amount of flow through the body due to the local dif- 


ferential pressures. 
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Figure 31. Differential Pressure for Cp = 4.45 at T = 8.85 
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Figure 32. Differential Pressure for Cp = 4.45 at T = 11.29 
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Figure 33. Differential Pressure for Cp = 4.45 at T = 12.92 
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Figure 34. Differential Pressure for Cp = 4.45 at T = 14.54 





Figure 35. Differential Pressure for Cp = 4.45 at T = 16.98 
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Figure 36. Differential Pressure History Comparison 
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IV. CONCLUSIONS 


The results presented herein warranted the following 


conclusions. 


lu. 


The effect of porosity is a function of the differential pressure 
across the porous material, and is dictated only by the physical 
characteristics of the material. 


. For a given canopy, the normalized differential pressure (the 


pressure coefficient) is independent of the ambient velocity, up 
to a velocity of about 600 miles per hour. 


. The effect of porosity cannot be adequately accounted for by let- 


ting through the camber a constant uniform velocity, as done by 
Inoue. The effect of the variation of differential pressure with 
time and along the camber must be included. 


. A novel method has been devised which combines the uniform 


porous flow concept with the variability of the flow rate with dif- 
ferential pressure through the use of discrete vortices embedded 
in the camber. 


. The results have shown that, in general, the effect of porosity is 


to displace the wake vortices further downstream and to retard 
their return to the camber during the period of deceleration. 


. Even the most porous canopy material used in the construction of 


large parachutes is not sufficiently porous to prevent the collapse 
of the parachute, as evidenced by the field experiments and veri- 
fied by the present numerical model. 


. The results of calculations with a much larger porosity have 


shown that it will be possible to considerably delay or entirely 
eliminate the collapse phenomenon. In doing so, it may be 
advantageous to optimize the radial distribution of porosity along 
the camber in order to achieve as large a drag force as possible. 
One must, however, bear in mind the fact that one also needs to 
account for the flexibility of the canopy. This forms the basis of 
the continuing efforts toward the understanding of the behavior 
of aerodynamic deceleration devices. 
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